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Abstract 

Lattice gauge theory's discretization of spacetime suffers from a drawback in that Lorentz co- 
variance is lost because the axes of the lattice create preferred directions in spacetime. Smaller 
and smaller lattice spacings decrease the effect but fail to eliminate it completely. It has been 
argued recently that detecting such a set of preferred directions or similar constraints would in- 
dicate whether the universe itself has an underlying lattice, i.e. the digital universe hypothesis. 
In this paper, I demonstrate a technique for accomplishing lattice gauge theory simulations while 
■ maintaining exact Lorentz covariance by replacing the lattice with a lattice graph such that the 

metric is defined as a discrete, Lorentz covariant matrix potential over the graph rather than a 
metric over an underlying manifold. This technique eliminates the symmetry violation of standard 
lattice gauge theory and suggests that, even in a digital universe, Lorentz covariance can still hold. 

PACS numbers: ILlS.Ha, ILSO.Cp, 12.38. Gc 
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Standard lattice gauge theory was first proposed by Wilson in his seminal paper {g] on the 
confinement of quarks. In his computational paradigm, the spacetime manifold is discretized 
by a hypercubic lattice. Actions for theories such as QCD are computed by summing over the 
plaquettes (closed loops formed by four edges) of the lattice via the discrete Yang-Mills action 
for the plaquette. The form of the action as a product of symmetry group members maintains 
the symmetry group of the theory, e.g., SU(3). The symmetry group of spacetime itself, 
however, is not maintained. In other words, rotations, boosts, and translations of the lattice 
cause the potentials and sources and hence the action on the lattice to change. If the universe 
has Poincare symmetry (or de Sitter symmetry), then lattice gauge theory only approximates 
this symmetry for sufficiently small lattice spacings. Recently, it has been suggested that 
the universe may have an underlying lattice, the so-called "computer simulation" universe 
a where real phenomena a.e the .e.ult of =oph.ticated lattice gauge co,nputat.o,., and 
that this could be detected by violations of the Greisen-Zatsepin-Kuzmin (GZK) limit [2[, 
resulting in a breakdown of special relativity at high cosmic ray energies. The authors, 
however, admit that computational paradigms "beyond our comprehension" could be used 
to eliminate the violation they are looking for. Observations that constrain violations of 
Lorentz invariance include astrophysical jSj and neutron electic dipole measurements [1] but 
these are at much lower energies than the GZK. In this paper, I show that lattice gauge 
simulations need not violate Lorentz symmetry and that the computational method to do so 
is well within the realm of human comprehension. Therefore, even in a computer simulation 
universe, Lorentz symmetry is not necessarily violated. 

Lattice gauge theory violates Lorentz invariance because the lattice is embedded in a 
mathematical abstraction, the manifold, which represents "real" spacetime Q. The manifold 
is continuous and differentiable and potential and source fields used in physics computations 
are defined as differentiable functions on it. Lattice gauge computations represent the lattice 
explicitly as sequences of vertices connected by edges while the existence of the manifold is 
implicitly understood. Because the computational physicist chooses the lattice's orientation 
within the manifold, a non-physical "preferred" coordinate system is imposed on spacetime. 
For computational purposes, as long as the lattice spacings are small, the coordinate systems 
should not affect the computations. This is not guaranteed in all cases, however, because, 
as in the growth of crystals, the cubic structure can affect the "macroscopic" scale under 
the right circumstances. As lattice gauge computations become more complex with greater 



2 



computing power, e.g. Petaflops (such as with DOE's Jaguar supercomputer and the up- 
coming Titan) and Exaflops hkely to be achieved in the next decade, problems may arise 
because of the underlying cubic crystalline structure of the lattice. (DeGrand, et. al j4| is a 
good reference for lattice gauge methods for QCD.) 

To solve this problem, an alternative to the standard lattice approach is to define the 
lattice with no underlying manifold so that it is just a graph. To be clear, when there is no 
underlying manifold, the edges and vertices have no position (or length) whatsoever, and 
rotations, boosts, and translations of the lattice itself have no meaning. Potential and source 
functions become discrete functions on the vertices, which relate to one another via the edges 
of the graph. Hence, the cubic crystalline structure of the former preferred coordinate system 
no longer exists because there is no coordinate system. 

Removing the manifold creates several new problems. There is now no coordinate system 
and hence no way to measure distances between event points that occur at the vertices, 
which is essential to quantum field theory. There is also no differentiability, which means 
that the continuum limit Lagrangians that contain derivatives cannot be obtained from the 
lattice as it now stands. Intuitively, it seems like a manifold is essential to defining a theory 
that describes quantum behavior, but, as 1 show in this paper, it is not. Both of these 
problems are solved, instead, by defining a gauge theory on the lattice graph such that 
Poincare or de Sitter covariance is conserved, and the continuum limit recovers the standard 
Yang-Mills action. The gauge theory ensures that the lattice reflects all the symmetries of 
a fiat manifold without the manifold. No curved spacetime manifolds are considered. 

First, an one dimensional example distinguishes between the standard approach and the 
approach I take below. Consider the real number line on which is defined a metric. For 
an arbitrary continuous, differentiable function defined on the real numbers, / : R — )■ M, 
a translation invariant functional S on /, S[f] = J dxL[f{x)] can by approximated by 
evaluating it over a lattice Co = {. . . , — 2e, — e, 0, e, 2e, . . . }, S[f] ^ J2h=-oo <^-^[/(^^)] = ^e[f]- 
If we translate the lattice: £o = {• • • , — 2e, — e, 0, e, 2e, . . . } — )■ £5 = {. . . , —2e+6, —e+6, 6, e+ 
(5, 2e + (5, . . . } by an amount 6, the approximation to the functional is, for arbitrary / and 
S, slightly different than the untranslated approximation: S^^s[f] = S^[f] + a for some 
small difference a. Therefore, the functional approximated on the lattice is not translation 
invariant, unlike the functional on the complete real line. In the continuum limit as the lattice 
spacing vanishes, however, translational invariance is recovered: S^^o^slf] = S^^o[f] = S[f]. 
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In higher dimensions, we also have boosts and rotations in addition to translations, but 
the principle is the same: discretizing over a lattice within a spacetime manifold eliminates 
the natural covariance of the manifold. 

The translation invariance can be added by redefining the lattice as a graph. The 
lack of translation invariance started because we assumed the lattice was embedded in 
the real number line. Suppose we eliminate the real number line and define, instead, 
an infinite graph where each vertex has exactly two neighbors. Label the vertices V — 
{. . . , y_2, y-i,yo,yi,y2, ■ ■ ■}■ Define a function : V ^ M. on each vertex. This function 
is a distance function and gives the distance between yk_i and yk- Let /fc : y — > R be an 
arbitrary function defined on the lattice. The approximation to the functional above is given 
by Sglf] — Yl'k'=-Qo dk-^ifk] ^-iid is translation invariant because gk is translation invariant. 
Let gk^e and fk = f{ke) and Sg[f] = Se[f]. If fk satisfies 



9fc-*-0 gk 

finite such that it is differentiable in the limit, then Sg[f] — )■ S[f] as e ^ 0. Although this 
example is too simple in one dimension to demonstrate the gauge covariance of g^ which will 
become evident in higher dimensions, it illustrates the main feature of the approach, which 
is to eliminate the metric on the underlying manifold in favor of a distance function that is 
tied to the lattice. 

Extending this approach to higher dimensions is much trickier than simply defining a 
metric function gj^ that gives the distance between every pair of vertices yj and yk because 
it does not provide a way to represent coordinate transformations. In one dimension, there 
is no coordinate freedom in the metric. The distance function gk, which, in our special case, 
was simply e, is a scalar field and necessarily invariant. In higher dimensions, one gjk and 
another g'jk may be equivalent under a coordinate transformation, e.g. rotation, because 
they covary. In the standard, lattice-in-a-manifold approach, the manifold takes care of 
that covariance via transformation rules for differentiable manifolds. On a lattice with no 
manifold, different rules apply which are, nevertheless, equivalent in the continuum limit. 
These are the rules of gauge covariance: 

(A note on subscripts: The beginning letters of the latin alphabet a,b,c, . . . ,h represent 
4- vectors. These letters range over 0, 1, 2, 3. I will use y4 to refer to the additional, 5th 
dimension that completes de Sitter group vectors.) 
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First, we must define the gauge theory as in the 1-D example but now for a 3+1-D 
universe. Let time be imaginary t — )■ so that spacetime is Euchdean and boosts become 
rotations, and define the coordinate system as de Sitter. (This is for simphcity. The following 
theory is achieved in the same way for the Poincare symmetry group.) On a four dimensional, 
infinite lattice graph such that every vertex has eight neighbors, label each vertex with five 
numbers yo,yi,y2,y3,y4- These labels become the coordinate system once the metric is 
defined. The fifth component, y^, is constant under the Poincare group, with 7/4 = 1 for 
all vertices, but can change under the de Sitter group. For three vertices, y neighbor to z 
neighbor to w, let f/yw(z) be a sma// rotation matrix with e > the small parameter. The 
matrix Uy^ = exp{eAa) can be split into two real valued parts: Gab and Habc as a sum of 
generators: 

Aa = GabVb + ^HabcMbc, (l) 

where Vb are generators for rotations in the yb-y4, plane and Mbc = —M^b are generators for 
S0(4) for rotations in the yb-yc plane where a,b,c = 0,1,2,3. Each of the ten generator 
matrices, Vb and Mbc, is a 5 by 5 matrix. 

A note on terminology: These matrices generate rotations in all 10 planes of a five 
dimensional space, but the VJ, generators can be considered "translations" for sufficiently 
small rotations, where the de Sitter group approximates the Poincare group. Therefore, 
in the following, I refer to the transformations under Vb generators as "translations" with 
the understanding that, in the case of the de Sitter group, these are actually rotations 
with respect to a fifth coordinate axis. I refer to Gab as the metric potential and Habc as 
the torsion potential because of their roles in determining distances and torsions (twists in 
parallel transport of vectors). In general, torsion will be zero. 

Now, we need to extract a metric from our "spin-1" matrix potential, Aa. To linear order 
in e the rotational matrix splits into component form, 

Uy^ = I + eAa = I + e{GabVb + ^HabcMbc). 

Using this transformation, we can define the coordinate system in terms of the transfor- 
mation potentials. Let y = {yo,yi,y2,y3,y4) and w = {wq, Wi, W2, W3, W4) , then 

w = [/ + e{GabVb + ^HabcMbc)]y. (2) 

The vector y experiences, on its way to becoming w, a transformation (translation or rota- 
tion) by passing through z. The formulas ti VaVb = Sab and ti MabM^d = Sac^bc 3] eliminate 
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The first term in the right hand side of [3] represents the torsion or twisting of a vector as 
it travels from y through z to end at w and is included for the sake of completeness but 
will be set to zero later. The second term on the right hand side represents the amount of 
translation. (Note: if Gab{'^) = then there is no translation and the two lattice points are 
in the same location.) 

To see how the above potentials describe a coordinate system without reference to a 
manifold, consider an arbitrary matrix potential component Bi{ya) as an example. Label 
the eight edges attached to vertex ya in pairs as 1, — 1, 2, — 2, 3, — 3, and 4,-4. Then if 
Bi{ya) > 0, its direction corresponds to edge 1. If -Bi(?/a) < 0, we have edge —1. The 
edge's direction in coordinate space is defined by [3l Likewise, the position ya relative to its 
neighbors is defined byO 

Now, we need to obtain a standard Euclidean coordinate system. Given that spacetime 
is flat (gravity is insignificant), there exists a gauge transformation Gab G'^f^ such that the 
metric potential is Euclidean: = 6ab and if^^^ = 0, where, for the 5*0(5) transformation 
0{x) such that O^^ = O^, the gauge transformation of the matrix potential in[T]is 



This equidistant Cartesian coordinate system corresponds to a lattice where each vertex 
is a distance e/2 from each of its eight neighbors. Note: the distance between two non- 
neighboring event points is undefined and is not necessary to lattice gauge simulations. 
Instead, relationships between non-neighboring points are given indirectly by probabilistic 
correlations, which the lattice gauge simulation samples. In the classical limit h 0, these 
correlations lead to the "macroscopic" notion of distance, but this is beyond the scope of 
lattice gauge theory. 

Now, we define our lattice gauge theory. The metric and torsion potentials are defined on 
the vertices between the event points. There are three kinds of vertices: event points where 
sources are defined, transition or transformation points where gauge potentials are defined, 
and action points where actions are defined (Fig. [1]). This definition is equivalent to that 
given in j?!. 



A'^ = OAaO^ - {daO)0^. 



(4) 
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FIG. 1: The plaquette consists of three kinds of points: event points, transition points, and an 
action point. This ensures that every quantity relevant to the lattice gauge computation is defined 
on a vertex and avoids having the action, in particular, appearing to be defined beyond the graph. 

In the Wilsonian derivation of any Yang-Mills theory, e.g. QCD, the metric appears from 
summing over plaquettes 

aa. 

Given an action as a sum of products of transformation 

matrices, 

S,=tT Yl U,U,UkUi, (5) 

plaquettes 

where each plaquette has sides A;, / in a closed loop and f/j, Uj, Uk, Ui G SU{N) 5*0(5). 
Each matrix has an SU(N) submatrix in the upper left corner and a constant S0(5) sub- 
matrix in the lower right corner. Eqn. [5] is invariant under a change of coordinate gauge 
Va X^aUb + i^a bccausc, by [31 the potential values defined on the vertices are fixed irre- 
spective of how the coordinate system changes. 

Since the method is restricted to straight-line coordinates such as Cartesian coordinates 
and global, i.e. constant, S0(5) transformations, the usual gauge transformation reduces to 
an ordinary rotation transformation: 

A', = OA^O'^, (6) 
such that, in real time, ordinary Lorentz transformations transform the metric: 

and torsion remains zero: Habc = 0. Thus, in equidistant Cartesian coordinates. Gab = Vab- 
Because the S0(5) part of each U matrix in [5] is constant, it also drops out of the final 
action leading to a continuum limit identical to standard Yang-Mills theory. By [3l the 
metric potential Gab appears everywhere the manifold metric usually appears. 
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Because Gab is the same as the standard metric for a flat spacetime, lattice gauge simu- 
lations with the modifled method will be identical to those with the standard method for a 
given orientation of the lattice. If the lattices are transformed by a Lorentz transformation, 
however, a difference between the two will arise. The difference between the background 
manifold and no background manifold methods is illustrated in Figures [2] and [3] for a transi- 
tion between two frames which are rotated with respect to one another by a transformation 

from a frame K to K'. 




FIG. 2: The plaquette with a background manifold transitioning from frame K to K'. 
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FIG. 3: The plaquette without a background manifold transitioning from frame K to K'. 

This paper has shown a subtle difference between a lattice embedded in a manifold and a 
lattice with no manifold. In the absence of gravity, there is no difference between the two in 
the continuum limit (assuming differentiability holds on all functions over the lattice), but, 
before the continuum limit is taken, only one is Lorentz covariant. This shows that even in 
a so-called digital universe of discrete lattice points, observables can be absolutely conserved 
under continuous Lorentz transformations. Further work would attempt to show adverse 
effects from Lorentz symmetry violations in lattice gauge simulations and also extend to 
curved spacetimes. 
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